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Abstract 

The evolution of matter produced in nucleus collisions is investigated in the frame 
work of the simplified fluid equation, where the compressible bag model is used 
as the equation of state for hadrons. It is shown that the quark-gluon phase is 
easily achieved in RHIC experiments, and that the typical structure of the phase 
transition, however, may be difficult to observe in simple analyses of (pt) data. 



*e-mail: kagiyamaOsci .kagoshima-u .ac.jp 
^e-mail: minaka@edu . kagoshima-u .ac.jp 
■fe-mail: nakamuraOsci .kagoshima-u .ac.jp 



1 Introduction 



There is long expectation to find the signal of quark-gluon plasma formation in the high 
energy experiments of nucleus-nucleus collision. Indeed, not a few of the experimental results 
suggest the existence of the quark phase but it seems that the incident energy is a little 
lower to observe the clean signals even at CERN SPS. The recent experiment at BNL RHIC 
may give us cleaner and richer signals. 

The purpose of the present paper is to give a manageable model for evolution of hadronic 
matter and quark-gluon plasma, produced in nucleus collisions, and to give some numerical 
predictions. There are two main ingredients in the model; the evolution equation of the 
matter and the equation of state for hadrons. 

In the evolution equation, assumed is cylindrical symmetry with respect to the collision 
axis, and used are the averaged quantities over the transverse plane. The averaged evolution 
equation is already reported in || ||, and its benefits are in that the parameterization 
of the unknown initial-conditions becomes very simple and in that the concerned physical 
quantities, such as size and transverse velocity of the matter, are directly obtained, in contrast 
to the straightforward calculation of the relativistic fluid equation M. 

As for the equation of state for hadrons, the compressible bag model is used, 
which has nice features as explained in the following. The simplest point-like pion gas 
model is often used but is lack of degree of freedom, and then it cannot reproduce the 
observed multiplicity of hadrons as long as natural parameters are employed. On the other 
hand, if the higher resonances are fully introduced in the point-like model, there appears 
the limiting temperature and the deconfinement transition to quark phase never takes place. 
The compressible bag model gives a solution for this dilemma, which gives sufficient degree of 
freedom and sufficiently high temperature because the mass of the hadrons becomes heavier 
as they are compressed by surrounding pressure. In this model the phase diagram with 
baryons is also calculable explicitly for whole region of parameters, and it coincides with the 
expected shape 0. 

In Fig. 1, shown are the three steps in the nucleus collisions: At the first step, gluons 
are exchanged between the nucleons, and the both nucleus are shifted to higher color-charge 
state. The initial energy density is essentially given at this stage. In the second step, the 
color flux is expanded in the longitudinal direction, and at proper time tq ~ 1 fm, the 
flux materializes into hadrons, quark-gluons or their mixed phase depending on the energy 
density. In the third step, the matter expands both in longitudinal and transverse directions, 
being cooled down to the freeze-out temperature, and then breaks into individual hadrons. 
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Fig. 1: (1) Step 1. Gluon exchange between nuclei. (2) Step 2. Longitudinal expansion of 
the flux tube. (3) Step 3. Cylindrical expansion of the hadronic or quark matter. 

Since the present model treats the third step, in order to show explicit results of the 
model, we must take some assumptions for the initial energy density of the produced fluid. 
Therefore the details of the results depends on the assumption, but it is shown that the quark- 
gluon phase is achieved in RHIC experiments and that the signal of the phase transition, 
however, is difficult to observe in simple analyses of (pt) data. The latter is mainly because 
of the richness of freedom in hadron phase: The richness pushes up the entropy density in 
hadron phase, and its difference between hadron and quark phases become much smaller 
than that calculated in simple pion gas model. Therefore the period of mixed phase becomes 
also smaller, and the typical structure of the first order phase transition is smeared out. 

In the next section, the equation of state in the compressible bag model is briefly reviewed 
on the basis of the results of 0. In §3, evolution equations are presented for perfect fluid 
of hadron or quark-gluon or their mixture and it is revealed that the evolution equations 
have an expanding solution as well as a scaling solution. And relations to observables are 
also given. In §4, the results of numerical calculations are given and discussed on the basis 
of the expanding solution for iViV-collisions and v4v4-collisions. Final section is devoted to 
conclusions and discussion. 

2 Equation of states 

In this section, we briefly summarize the results of [|] [0]. Let us suppose a mixed gas of 
n- kinds of meson enclosed in volume V at finite temperature T. In the compressible bag 
model, the free energy function F of the gas is given by 
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where N{ is a number of z-th kind of meson and Vi is its volume and Mi its mass. The constant 
b is a volume exclusion efficiency parameter. The function Ff is a free energy function of 
free (point-like) boson gas. As for mass function Mj(fj), we assume that of MIT bag model: 



Mi(vi) = A iVl 1/3 + Bvi, (2.3) 



where B is bag constant. 



Under the approximation that the average of inverse Lorentz factor of z-th kind of meson 
equals to unity ((7j rl ) ~ 1) in the rest frame of the system, basic requirement of the com- 
pressible bag model that dF/dvi = and the requirement that the the chemical potential of 
meson should vanish, /ij = dF/dNi = 0, determine the pressure p of the system as a function 
of T by the following equation: 

d 3 k 

W) 



5> T J ^^{l-eM-(Ei-ri( P ,rni))/T}}, (2.4) 



p = - 

E i = v /fc 2 + M i (p,m i ) 2 , 

M l (p,m t ) = mjl + ^j + , (2.5) 

biriij) ( bp \ 

Hiip,™*) = ~bvip= — T^"( 1 + ^"j ' ( 2 - 6 ) 

where gi is a degeneracy factor of the z-th kind of meson and mj = 4(Aj/3) 3//4 S 1//4 is its mass 
in the vacuum. 

Now we take the limit n —>■ oo and replace the sum with the integral 

n . 

Y.9i=> J dmr{m), (2.7) 
i=i 

where r(m) is a level density function. Then ( [2.4|) is written as 

d 3 k 



P 



- J dm r H T / ln I 1 - e- (s ^' )/T ] , E = Vk 2 + M 2 , (2.8) 
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Thus we can calculate p as functions of T from ( |2.8|) when an appropriate level density 
function is given. An entropy density s is given by 

dp C . 
S = W = D< (2 ' n) 



3 



V r , , n r d 3 k (E-u')/T , 

C = - + dmr(m) - — - £ 2.12 

T J V J J (27r) 3 exp[(E-/i / )/T] -1' V 1 

D = 1+ / dmr{m) f f^—, ^7 ^ , (2.13) 
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(2.14) 



An energy density e is calculated by the following thermodynamic relation 

e = Ts-p. (2.15) 
A number density p for the hadron with vacuum mass m is given by M 



P' 



1 + 6/ dmr(m)p'v ' 



(2.16) 



/ / d 3 k 1 
P = 7 (27r) 3 exp[(£ -//)/T] - 1' (2 ' 1§) 

Precisely speaking, the above p is not the number density but one should multiply it by r(m) 
and integrate over some mass range, say [m — Am, m + Am], to get the number density of 
the hadron with vacuum mass around m. In practice, however, we use the above p simply 
multiplied the degeneracy factor of the hadron, for example 3 for pion, instead of integrating, 
as the number density. This prescription is justified when narrow width approximation is 
applicable, which we will assume in this paper where only pion density is considered in 
numerical calculations. 

As for the level density, we adopt the following form 

r(m) = 35(m — m^) + a9(m — 2m 7T )m~ d exp(m/Tu) , (2-19) 
and whose parameters are 

a = 25.462^, d = 3, T M = (9(LB/167r 2 ) 1/4 , (2.20) 
according to the Abelian bag model by Kapusta M. 



3 A simple fluid model 



4 



3.1 Evolution equations 



In this subsection, we present evolution equations of perfect fluid of hadron or quark-gluon 
or their mixture. In the previous paper f3|, hadrons had been assumed to be massless, but in 
the present paper, hadrons are assumed to obey the equation of states given in the previous 
section. 

As shown in §3 of |§, transverse flow with cylindrical symmetry obeys the following 
equations: 

T = I 3 - 1 ) 

GST T 

f = -'-±2, (3-2) 
dr r 

where r is defined by (t 2 — z 2 ) 1 ^ 2 , and has the meaning of the proper time only on the collision 
axis (z-axis). s, e and p are quantities integrated over the transverse plane; 

s = J d 2 r sjt, e = J d 2 r [(e + p)^ 2 — p], p = J d 2 rp, (3.3) 

where jt is the transverse 7-factor of the fluid, and e, p and s are energy density, pressure 
and entropy density respectively. 

Following , let us introduce a radius R of the cylindrical fluid and assume the following 
kinematical relation; 

^ = v t , (3.4) 

where Vt is the transverse velocity and is related to 7* by the relation 74 = (1 — v 2 )~ 1 ^ 2 . And 
further let us take a crude approximation, 

s = 7nR 2 S7t, i = tcR 2 [(s + p)i 2 t - p], p = rcR 2 p, (3.5) 

and hereafter regard s, e and 7$ as the averages over the transverse plane, or as the values 
at the representative point of r. All the quantities, then, only depend on r. 

Integrating fl3.1|) and using (|3.5|) , we get the following equation 

ts = 7iR 2 S'-yt T = const.. (3.6) 
Differentiation of this equation with respect to r yields 

2R s 7 t 1 
= — + - + - + - 

R S J t T 



2v t s v t 1 



+ Z + T-^Vt + -, (3.7) 



where dot represents differentiation with respect to r. 

Before proceeding to derive an evolution equation, we should note the following relations 
obtained from the thermodynamic relations ( j2.11|) and ( 2.15| ), 
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(3.9) 



«t = 0, 



(3.10) 



or 



T 



Vt T + T 



2(1 - ^ 2 )p 
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(3.11) 



vl R(e + p) 

The equation ( |3.10| ) has a scaling solution R = R , s = s t /t, where R and s are initial 



values at r = tq. As for (|3.11| ), it together with 



v; 
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f 
T 
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is recasted as 

2 P {i-vfy 



1 — at;/ 
s 



2 ih 



e + p 

P(l-«?) j , 1 
e + p J r 



(3.12) 

(3.13) 
(3.14) 



Tds / dT 

The equations (|3.4| ), (|3.12|) and (|3.13| ) are a closed set of evolution equations of fluid, which 
is also applicable for the fluid in quark-gluon phase as well as hadron phase. These equations 
have an expanding solution if an initial value v t (ro) is non-negative. In the next section, we 
will numerically analyze on the basis of this expanding solution. 

For quark-gluon phase, equations of states are given by 



P 



K Q T* - B, 



4:KqT , e = 3K Q T* + B 



(3.15) 
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with 



(7N* 8 \ 2 

Kn = 1 H ) 7T 

Q V 60 45 



(3.16) 



where Nf is a number of quark flavor. The evolution equations ( |3.12|) and ( p,13|) for quark- 
gluon fluid read 



f 
T 



l^g U 3 + vj - 3(1 -vf)B/K Q T* 
3 - v? \ r 2R 
1 fv t (5-v*) , ll 



2i? 



(3.17) 
(3.18) 



For mixed phase, equations of states are given by 

P = Pc = Pq(T c ) =p H (T c ), 
s = xs Q (T c ) + (1 - x)s H {T c 
e = xe Q {T c ) + {l-x)e H {T c 



(3.19) 
(3.20) 
(3.21) 



where T c is critical temperature and the subscript H (Q) denotes the quantity for hadron 
(quark-gluon), and x is a fraction of quark-gluon. Corresponding equations to (|3.12|) and 
( [3.13j ) are given by 



v t 



x 



2p (1 - n 
: + p R 
s 



2\2 



s Q (T c ) - s H {T c ) s Q (T c ) - s H {T c 



2v t I p(l-v 2 t ) 
R \ e + p 



1 

+ - 

T 



(3.22) 
(3.23) 



Thus we have obtained evolution equations for all phases. Those are (|3.4|) , ( p,12|) and ( p,13|) 
for hadron phase, (fOj), (|3.17| ) and Q3.181 ) for quark-gluon phase and (|3~^ ), ( p.22j ) and ( p.23|) 
for mixed phase. Throughout the evolution, ts is conserved, as is seen in ( |3.6| ), so that it 
can be used as a check in solving the evolution equations numerically. 



3.2 Relation to observables 

Now that we have the evolution equations, we can solve them and get T (x for mixed phase), 
R and % as functions of r if some initial values T , (a;o for mixed phase) and i?o are given 
and the initial condition 7t(ro) = 1 is assumed. The final condition of the evolution is set 
by the freeze-out of hadrons from hadronic fluid. If the freeze-out temperature Tf is given, 
we can get 77 from the solution, 

T f = T{jf) , (3.24) 
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where the subscript / means freeze-out. Then we get the values of the radius R and 74 at 
the hadronization; 

R f = R(r f ), ltf = lt (r f ). (3.25) 
These quantities are related to experimental observables. 

(i) Charged multiplicity 

The produced hadron density in rapidity space is given by 

= Tf^Rjpntf, (3.26) 

where p/ is the number density at r = 77. If we assume that all the produced hadrons are 
pions and put N c ^ = (2/3)iV, we get the charged particle density, 

dN ch 2 „ 

-— = -TfirR f p v fytf, (3.27) 
where p n denotes effective pion number density defined by 

Ptt = Pn + NmPm, (3.28) 

Pn = = = - , /* rT7 - , (3.29) 
1 + J amT{m)p'v 

Pm = r / M 1 v ; , (3.30) 
1 + J dmr\m)pv 

1 f d?k 1 

P?r = 3 i (2 7 r)3exp[(K- / ,0/T]-l' (3 - 31) 



Pm = / rfmr(m) /" — (3.32) 

J2 m7r j (27r) d exp[(i? — p )/T\ — 1 



In (|3.28| ), the first term corresponds to pole term and the second term corresponds to con- 



tinuum part of r(m) in ( [2.19| ). The quantity Nm is an average multiplicity to pions from 
the continuum of meson resonances. 

(ii) Transverse momentum and transverse energy 

The average transverse momentum of pion in the center of mass system is approximately 
given by 

<J?t) = >J(mfh?-ml, m* t = yjpf + ml, (3.33) 

where pi (^1*) is a transverse momentum (mass) of pion in the rest frame of the fluid. The 
average transverse mass of pion in the rest frame of the fluid is calculated as an average 



of the one for directly produced pion and the one for the contribution from pions as decay 
products of continuum of meson resonances; 



. ,2 , . P*f (mf M) + N M pM f (mf (tt c ) ) 



«> = rJVt —--^^ (3.34) 

At/ + ^MpMf 

The quantity (m^ 2 (7r c )) is estimated on the basis of random walk approximation as 

K 2 (7r c )> = (mf(M))/N M . (3.35) 
Then it follows from ( ggg ), (|3~30D , (ggg) and ( g3| ) 



} = *, + jw„, • < 3 - 36 > 

d 3 fc A; 2 + M 2 



K 2 W) = (3.37) 

(27r) 3 exp[(^- yU ')/^]-l 

d 3 fc A; 2 + M 2 



dmrim) 



m?(M)) = 1*= 7 (yexp^-ziO/ri-l 

* v y/ r 00 , , , /• d 3 fc l 



/2m. v y 7 (2tt) 3 exp[(E 7r - /i')/T] - 1 
where fc t denotes a transverse momentum variable. 

Transverse energy per unit rapidity interval is simply given by 

dE t r. — — dN 

V = < 339 > 

(hi) Rapidity interval and total multiplicity 

Since the evolution of fluid does not take place yet at r = r , total energy of the fluid 
E tot in the center of mass system is calculated as 

/?7max 
t e itR cosh rjdr] = 2r e nR sinh 77 max . (3.40) 
"?7max 

Then the rapidity interval of the fluid 2r] max is determined when E tot and initial temperature 
T (x for mixed-phase start) are given, that is, 

2r/ max = 2 sinh- 1 ( - Etot ) . (3.41) 

Total multiplicity (N cil ) is approximately given by 

(N ch ) » 2 W ^. (3.42) 
ay 



9 



4 Numerical calculations 



4.1 AW-collision 



In order to start numerical calculations, we have to fix bag constant B, volume exclusion 
parameter b, the number of quark flavors Nf, average multiplicity of pions from continuum 
Nm, initial proper time tq, initial temperature To (xq for mixed-phase start), initial radius 
Ro and freeze-out temperature Tf. First of all, we fix as 

N f = 2. (4.1) 

The initial radius Rq for iVW-collision is set to the experimental value of proton charge radius 



R (NN) = 0.82 fm. 
The parameter b is determined by the relation 

bv N (p = 0) 



^ = ^(0.82fm) 3 , 
AB 3 v ' ' 



(4.2) 



(4.3) 



when the bag constant is given. We adopt the following values 

B 1/A = 0.25 GeV, 6 = 5.00. (4.4) 
Then critical values at the deconfinement transition are given by 

T c = 0.196 GeV, e H (T c ) = 0.874 GeV/fm 3 , e Q (T c ) = 2.85 GeV/fm 3 . (4.5) 

For the purpose of determining Tf and tq, we use the fit to experimental data of iViV- 



collisions as inputs [10 



Tf 

its 



dN. 



ch 



dy 



0.135 ±0.003 GeV, 

1 + (0.010 ± 0.003) log(E cm /E 1 )6(E cm - 
Ex = 12 ± 8 GeV, 

(1.31 ± 0.02) + (0.20 ± 0.03) \og(E cm /E 2 
+ (0.062 ±0.007) \og\E cm /E 2 ), 
E 2 = 12 GeV, 



(4.6) 
(4.7) 



where E cm is the initial energy in the center of mass system. We choose the central value 
0.135 GeV of (|4.6|) as Tf. Since the fluid immediately freeze-out at E cm = 12 GeV as seen 
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Fig. 2: (p t ) is shown as a function of dN^/dy. The parameter Nm is searched for the range 
2.5 ~ 3.0 with B l / A = 0.25 GeV. The experimental data are taken from (Tl| . 



from ( |4.7| ), the value of r is fixed to reproduce the value dN ch /dy = 1.31 at E cm = 12 GeV 



as seen from 

^ = -ToTTit^ = 1.31. (4.9) 

dy 3 

Now varying T (or x ), we can calculate (p t ) as a function of dN ch /dy with giving the 
value of Nm- The results are shown in Fig. 2. The value of Nm is searched for the range 
2.5 ~ 3.0. The value of (p t ) is almost determined by (ml 2 (ir c )) = (mj 2 (M)) /Nm since the 
value of N M pMf is almost an order of magnitude larger than that of p n f. (See ( |3.33| )^ (|3.35|) ) 
Therefore the calculated values of (pt) is sensitive to the value of Nm- Comparing with the 
experimental data [|TTJ we choose as 



N M = 2.7, r = 1.29fm. (4.10) 

The response to B is shown in Fig. 3. For small bag constant, (p t ) becomes large slightly 
fast as dN c h/dy becomes large. This is recognized as a consequence of a larger acceleration 
for smaller B as is seen from Q3.17| ). In Fig. 4, the initial temperature To is shown as a 
function of dN^/dy. For B 1 / 4 = 0.25 GeV the critical temperature is T c = 0.196 GeV and 
the region 4.2 < dN^/dy < 11 corresponds to the mixed-phase start. For a larger (smaller) 
value of B, T c becomes higher (lower) and the region of the mixed-phase start becomes 
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Fig. 3: (p t ) is shown as a function of dN^/dy with different B. The parameter Nm is set to 
2.7. For small B, (p t ) becomes large slightly fast as dN ch /dy becomes large. 




1 10 100 

dN ch /dy 

Fig. 4: T is shown as a function of dN^/dy. The parameter is set to 2.7. The critical 
temperature becomes lower and the region of mixed-phase start becomes narrower for smaller 
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Fig. 5: dN c h/dy is shown as a function of E cm . Solid line show a consequence of ( [4.1 1| ) and 
( |4.12| ) while dashed line show (|4.8|) . Both lines are coincide below Tevatron region 



wider (narrower). The value T c = 0.237 GeV for B lj/4: = 0.30 GeV is somewhat larger than 
the value suggested by lattice QCD calculations so that we have adopted _B 1//4 = 0.25 GeV 
in ( |4.4|) . Comparing with Fig. 3, we can see a slight structure of entrance to mixed phase 
while entrance to quark-gluon phase does not show a remarkable structure. 

An initial energy density Eq has a fundamental importance to prescribe an initial state 
of the fluid and determine its evolution in our picture as explained in §1. The values of Eq 
can be obtained as a function of E crn in such a way that the calculated values using ( ^.27[ ) 
reproduce the predicted values by (|4~8| ). However, ( [4.8| ) is only applicable for E cm < 1.8 TeV 
since ( |4.8|) was obtained as a result of fit to experimental data up to Tevatron region. In 
order to extrapolate above Tevatron region we assume the following form of Eq 

s (NN)=e 00 (^) 5 , (4.11) 



v E 2 , 

and determine Eqq and S to reproduce ( [4.8| ) below Tevatron region. The result is given by 

£ 00 = 0.184 GeV/fm 3 , 5 = 0.294. (4.12) 

Using this result dN c ^/dy is shown as a function of E cm in Fig. 5 where solid line show a 
consequence of ( [4. 11[ ) and ( f4.12p while dashed line show ( [4.81 ). We can see both lines are 
coincide below Tevatron region. 
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Fig. 6: dE t /dy is shown as a function of A. The parameter (3 is searched for the range 
of 0.40 ~ 0.45. E cm is set to 20GeV/A. Other parameters are set as B 1/4 = 0.25 GeV, 
Nm = 2.7. Experimental data are taken from [O. 



4.2 AA-collision 

First we set the initial radius Rq for AA-collision as 

Ro{AA) = 1.2A 1/3 fm. (4.13) 

In order to incorporate AA collisions into numerical calculations, we need a relation to 
combine Eq(AA) with e [NN). For this purpose we assume that 

£ Q (AA) = A p e (NN). (4.14) 

In the above, both sides should be compared at equal E cm which means energy per nucleon 
in the center of mass system for AA-collision. 

In Fig. 6 the average transverse energy per unit rapidity interval dE t /dy is shown as 
a function of A for E cm = 20GeV/A and the values of /3 = 0.40 ~ 0.45. Comparing 
with experimental data ||12| , where the initial energies are E cm = 19.4GeV/A (i?Lab — 
200 GeV/A) and E cm = 17.3 GeV/ A (E Lah = 158 GeV/A), we choose 

P = 0.43. (4.15) 

Having fixed all parameters we can now calculate various observables as functions of E cm 
and A. 
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Fig. 7: dE t /dy is shown as a function of A for various E cm . The parameters are set as 
B 1 ' 4 = 0.25 GeV, N M = 2.7, (3 = 0.43. Experimental data are taken from pi . 
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Fig. 8: To is shown as a function of dN^/dy. The parameters are set as B 1 / 4 = 0.25 GeV, 
Nm = 2.7, (3 = 0.43. The thin lines show the case of no phase transition for comparison. 
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Fig. 9: (p t ) is shown as a function of dN^/dy. The parameters are set as B 1 ^ = 0.25 GeV, 
Nm = 2.7, (3 = 0.43. The thin lines show the case of no phase transition for comparison. 



In Fig. 7 dEt/dy is shown as a function of A for E cm = 20 GeV/ A, 200 GeV/A, 2 TeV/A. 
dE t /dy grows almost linearly with A. In Fig. 8 the initial temperature T is shown as 
a function of dN ch /dy for A = 1,30,200. These results correspond to the experiments 
for various incident energies. For the case of A = 1 (Rq = Rq(NN) = 0.82 fm) 4.2 < 
dN c h/dy < 1.1 x 10 corresponds to mixed-phase start. For the case of A = 30 (Rq = 3.73 fm) 
8.6 x 10 < dN^/dy < 2.3 x 10 2 corresponds to mixed-phase start. As for the case of A = 200 
(Rq = 7.02 fm) dN c h/dy < 8.2 x 10 2 corresponds to mixed-phase start. The behaviour of T 
in the case of no phase transition is also shown by thin lines for comparison. 

In Fig. 9 the average transverse momentum (p t ) is shown as a function of dN^/dy for 
A — 1, 30, 200. These results also correspond to the experiments for various incident energies. 
The behaviour of (p t ) in the case of no phase transition is also shown by thin lines for 
comparison. It is clearly seen that the fluid is not accelerated so much in the mixed phase. 
The reason is that the acceleration in the mixed phase given by the right hand side of 
( p.22| ) has no 1/r-term and an extra factor (1 — vf) in l/i?-term so that it takes smaller 
values than the acceleration in hadron/quark-gluon phase given by ( |3.12 )/ (|3.17 ). This is 
a common feature to the case of massless hadron, which is supposed to be pion, discussed 
in the previous paper 0. The difference lies in the fact that the interval of mixed-phase 
start in the present case is narrower than the one in the case of a single massless hadron so 
that the structure is less prominent especially for A = 30, 200. This is because the degree 
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Fig. 10: (p t ) is shown as a function of dN^/dy with varying A — 1 ~ 240. The parameters 
are set as B l ' A = 0.25 GeV, N M = 2.7, (3 = 0.43. 



of freedom in the hadron phase is larger for the present case which takes into account the 
contribution from continuum of meson resonances. However it should be noticed that the 
entrance to quark-gluon phase can be seen in this figure plotted in the log scale of dN ch /dy 
in contrast to Fig. 3 plotted in the linear scale. 

In Fig. 10 (p t ) is shown as a function of dN ch /dy for E cm = 20 GeV, 200 GeV, 2 TeV with 
varying A from 1 to 240. To this variation the response of (p t ) is mild and does not become 
so large even for E cm = 2 TeV. In Fig. 11 the final radius Rf is shown as a function of 
dN ch /dy for E cm = 20 GeV, 200 GeV, 2 TeV with varying A from 1 to 240. The experimental 



data are taken from [13]. The calculated values does not differ so much from the observed 
values. It should be noticed that Rf is rather insensitive to the variation of E cm . As for the 
response to the variation of (3, Rf becomes larger for smaller f3. 

Finally in Fig. 12 the initial temperature is shown as a function of A for E cm = 20 GeV, 
200 GeV, 2 TeV. For E cm = 20 GeV the entrance to the mixed phase occurs at A = 27 and 
the entrance to the quark-gluon phase does not occur even at highest A. For E cm = 200 GeV 
the entrance to the mixed phase occurs at A = 6 and the entrance to the quark-gluon phase 
occurs at A = 87. As for E cm = 2 TeV the entrance to the mixed phase occurs at A = 2 and 
the entrance to the quark-gluon phase occurs at A = 18. Thus quark-gluon phase will be 
realized at BNL RHIC energies while it is not at CERN SPS energies. 
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Fig. 11: Rf is shown as a function of dN^/dy. The parameters are set as B 1 ^ = 0.25 GeV, 
Nm = 2.7, (3 = 0.43. The experimental data are taken from |L3 . 
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Fig. 12: To is shown as a function of A. The parameters are set as B 1 ^ = 0.25 GeV, 
N M = 2.7, = 0.43. 
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5 Conclusions and discussion 



We have shown that our simplified fluid model easily calculates observables in nucleon- 
nucleon and nucleus-nucleus collisions at high energy region. As for the equation of state for 
hadron phase, we have adopted the compressible bag model. This model can satisfy both 
requirements that one should take into account an infinite degrees of freedom of hadrons 
and that deconfinement transition should be take place, which usually contradict mutually 
in models other than compressible bag model. For an evolution of the fluid, we have assumed 
the cylindrical symmetry with respect to the collision axis and solved the averaged quantities 
over the transverse plane. 

In the previous paper |3j of the averaged evolution equation, we used the equation of 
state of massless pion for hadron phase. Both calculations give the result that agrees well in 
terms of qualitative analysis. However, the present approach produces (p n f + NMPMf)/p-nf 
times of pions in comparison with the previous one. This gives a reasonable value for To by 
the condition of ( |4.9|) . On the other hand, the ratio £q(T c )/£h-(T c ) becomes smaller because 
of the increase of hadronic degrees of freedom. 

The increase of the average transverse momentum is brought about by a transverse flow 
of the fluid. As is seen from (|3.12|) the pressure slope p/R in the transverse direction mainly 
causes the transverse flow. In the mixed phase the effect from the corresponding term 
became smaller as shown by equation QcS.lTj ). This gives a structure of an average transverse 
momentum as a function of dN^/dy. In the hadron and the quark-gluon phase the ratio of 
the pressure to the energy density increases rapidly with the increase of the energy density. 
The ratio in both phases reaches 1/3 in high energy density limit, so that the acceleration is 
generated similarly in both phases as shown in Fig. 9. In the case of the quark-gluon phase 
start the acceleration becomes weaker because of passing through the mixed phase. As for 
AA-collisions the acceleration becomes weaker with an increase of A because of a decrease 
of the pressure slope in the transverse direction. 

As shown in Figs. 8 and 9 dN c ^/ dy is not a reasonable variable for plotting the observables 
in AA-collisions except for the transverse source radius. In our model the equations of ( ^.12| ) 
and ( 4.14 ) generate the A dependence. We have to find a scaling variable in /1/1-collision. 

We have assumed the relation between the initial energy densities in NN- and AA- 
collisions by equation Q4.14J) . This power-law assumption is similar to the relations given 
by the dilute-gas and multiple-collision models ||. Finally we predict that the quark-gluon 
plasma can be observed for average multiplicity events in AA-collision with A > 87 at RHIC 
energies. The numerical calculation results are shown in Fig. 12. 
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